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ABSTRACT 



Propagation of electromagnetic energy through the atmosphere is a difficult task 
because of temperature fluctuations and index-of-refraction inhomogeneitics which de- 
grade the beam's coherence. Understanding this phenomena is of practical importance 
for optical systems. 

This thesis presents an analytical numerical technique which simulates the effects 
of atmospheric turbulence. The extended Huygens-Fresncl principle was used to simu- 
late wave propagation in a two-dimensional randomly varying medium, which is repres- 
ented by thin, filtered, Gaussian phase screens. The wave optics code implements both 
Fresnel and Fraunhofer propagation, by employing the fast Tourier transform (FIT) 
algorithm. The analytical spatial coherence length, p 0 , and normalized intensity vari- 
ance. G-jP, of the perturbed electric field, were examined. Results support the concept 
of intensity saturation for weak scattering cases, however, differences in the values of the 
theoretical and analytical spatial coherence lengths, occurred. 
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I. INTRODUCTION 



Atmospheric turbulence degrades the coherence of electromagnetic energy propa- 
gating through the atmosphere. The refractive-index inhomogeneities associated with 
the turbulent atmosphere induce phase and intensity perturbations across the wavefront. 
Understanding the propagation and scattering of optical waves in random media, is es- 
sential for atmospheric laser beam propagation and imaging systems. 

This thesis models the propagation through a random media by means of the ex- 
tended Huygens-Fresnel principle. A two-dimensional, thin, Gaussian phase screen re- 
presents the randomly varying medium. This wave optics propagation code employs the 
fast Fourier transform (FFT) algorithm in both Fresnel and Fraunhofer forms of prop- 
agation. The Fresnel region incorporates two forms of the diffraction integral, the 
transfer function form for "near field" distances, and the convolution form for "far field" 
regions. 

Previously, numerically intensive simulations of this type required large, super- 
computer systems. However, these computations were performed on a compact, 
desktop computer system. 

The model's accuracy was assessed by comparing computer values of the spatial 
coherence length, p 0 . with values obtained from the analytical mutual coherence func- 
tion (MCF). Additionally, the concept of intensity variance saturation was examined for 
a single phase screen in the Fraunhofer approximation. The normalized intensity vari- 
ance. o}//\ approaches a saturation value asymptotically for increasing values of the 
index-of-refraction structure parameter. Q . Theoretical calculations suggest saturation 
for multiple scattering, however, they say nothing about a single phase screen realiza- 
tion. This thesis provides results which support saturation effects for single scattering 
cases. 
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II. BACKGROUND 



A. STATISTICAL DESCRIPTION OF TURBULENCE 

1. Random Variables 

Maintaining the coherence of an electromagnetic wave propagating through the 
atmosphere, requires an understanding of the effects imposed on the wave by the tur- 
bulent medium. The fundamental characteristic of atmospheric turbulence is its ran- 
domness, which must be described statistically. In addition to statistical quantities, 
further assumptions must be made of atmospheric turbulence. These include the con- 
cepts of stationarity. homogeneity and isotropy. Stationarity implies that random 
processes are time independent. Homogeneity assumes invariance under a Galilean 
transformation of coordinates, while isotropic variables are invariant with respect to 
coordinate rotations. [Ref. 1] Mathematically, these two assumptions imply that the 
statistics at two points F, and r 2 depend only on the difference, r l2 = |F, — F,| . 

2. Local Homogeneity and Isotropy 

In general, atmospheric random variables do not obey the assumptions of 
stationarity. homogeneity, and isotropy. However. Tatarski [Ref. 2] introduces the 
concept of "local" homogeneous and isotropic random variables. This concept requires 
homogeneity and isotropy within a localized region of size L 0 , the outer scale length. 
Furthermore, the difficulties associated with nonstationary random variables, are re- 
moved by considering random fields with stationary first increments [Ref. 2]. Tatarski's 
method applies to a nonstationary random function whose mean varies slowly with time, 
by considering the difference of the function at two different locations. The slow func- 
tional changes do not affect the value of this difference. 

3. Structure Functions 

Tatarski introduces the structure function 

Djih - f,) = <W 2 ) -/(f,)] 2 >, (l) 

a tensor that is the difference between two quantities. Some important aspects of the 
structure function arc; that its general form is valid for any variable, and < > denotes 
an ensemble average taken over all possible point pairs r 2 . F,. Assuming homogeneity 
and isotropy, the vector dependence reduces to the magnitude r, 2 = |F : — F,| . and the 
structure function becomes 



